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We study the magnetism in the periodic Anderson model in the limit of large dimensions by 
mapping the lattice problem into an equivalent local impurity self-consistent model. Through a 
recently introduced algorithm based on the exact diagonalization of an effective cluster hamiltonian, 
we obtain solutions with and without magnetic order in the half-filled case. We find the exact 
AFM-PM phase boundary which is shown to be of 2 nd order and obeys ^jj- « const. We calculate 
the local staggered moments and the density of states to gain insights on the behavior of the AFM 
state as it evolves from itinerant to a local-moment magnetic regime 

PACS numbers: 71.27.+a, 75.20.Hr, 75.30.Kz 



I. INTRODUCTION 



The periodic Anderson model (PAM) consists of a band of d— electrons that hybridizes with localized /—electron 
states at each lattice site. The double occupation of the /—sites is disfavored by a local term that corresponds to the 
repulsive Coulomb interaction. 

This model hamiltonian is widely considered to be relevant for the description of a large class of strongly correlated 
systems, most notably, the so called heavy fermions and Kondo insulators. Examples of-tliese systems are, among 
I/*-) ' others, UPt 3 , CeAl 3 , the insulators CezBi^Ptz, CeNiSn, and possibly FeSi, and Smi^E'cl 

Despite much effort, our present theoretical knowledge about the variety of behaviors that solutions of this model 
may have, is still limited. Much objiiir current understanding of the model comes from approaches based in the 
consideration of a variational ansatzou, and also from, its formulation in limiting cases where the problem becomes 
y ■ more tractable, as for instance, the large N approachJjLl 

On the other hand, the steady increase of computational power, has recently made possible the numerical study of 
the model in 2-dimensional finite lattices by quantum Monte Carlo techniques. However, that numerical technique 
has same intrinsic limitations in regard of the parameter space that can be investigated, in particular the T — > 
£j , limit 

In this paper, we will consider the problem within the local imnurity self-consistent approximation (LISA). This 
. [ method was introduced in the present form by Georges and Kotliaia in the coptext of the Hubbard model, and later 
generalized for other models (including the PAM) by Georges, Kotliar and Si.E2l It can be thought of as a dynamical 
mean field theory for correlated electron systems. Its basic idea is to isolate a local site of the lattice model and embed 
it in an effective bath that contains the information of the response of the rest of the lattice. 

The LISA method becomes exact in the limit of large number of spatial dimensionslij Therefore, it also represents 
a limiting formulation of the model where simplifications occur (e.g. the self-energy becomes local), that allow us 
to obtain essentially exact results. This method has already bfen used in the periodic Anderson model for the 
investigation of the paramagnetic solution at finite temperature. E3 

It should be emphasized here, that due to the local character of the approach, it should be most suitable for 
the study of electronic systems where the relevant orbitals are very localized. Also, since the interactions can be 
considered in a non perturbative manner, it is most appropriate when the correlations play a crucial role. These 
features are clearly a characteristic of experimental systems such as heavy fermions and Kondo insulators where the 
PAM is considered to be most relevant. 

In this paper we focus on the study of solutions with magnetic order at T = 0. In particular, we will present 
the exact phase diagram. We will also study the nature of the different antiferromagnetic regimes that the solution 
displays. To this end we obtain the staggered magnetic moments and the densities of states in various regions of the 
model parameter space. I— ■ 

The question of the study of magnetically ordered solutions was previously considered by Doniachli-3 in the context 
of the Kondo lattice, and by EvanscJ and , more recently, by Sun et al£3 in the PAM. 
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We conclude with a brief discussion on the relevance that our results may have for the interpretation of experimental 
results. 



II. METHODOLOGY 

The formulation of the LISA method has been already discussed in detail elsewhere^ so here we shall limit ourselves 
to a brief presentation of the relevant expressions in particular we shall note how they are modified to consider solutions 
with magnetic order. 

We begin by writing the action of the lattice model in its Function Integral form 

5 = -E f f drdr' ^ ka {r)G^{r-r')^Ur') + f *r (n m (r) - \){n m ( T ) - \) (1) 

k,a • >0 ^° i •'° 

where the annihilation operator is defined as ipka = {fa, d^}, U is the local repulsion for double occupation of the 
jth y_ s ite, and k runs over the Brillouin zone. The inverse matrix propagator Gq 1 is explicitly 

where e° is the atomic energy of the /—site, /z is the chemical potential, and 14 is the hybridization matrix element 
between the / and the d— sites. Here, for simplicity, we take 14 = V, and consider the particle-hole symmetric case 
setting €j = /J, = 0. Since in the limit of d — > oo the self-energy becomes local, i.e. k— independent, we can write the 
inverse Green function as 



G a {kM = y v 



(3) 



To proceed further we consider the model on two sublattices A and B as is usually the case when dealing with 
antiferromagnetically ordered states, and demand the self-consistency condition that follows from the fact that the 
local Green function obeys 

1 f°° 

k J -°° 

We consider, for simplicity, the case of a Bethe lattice, which in the large dimensional limit has a semi-circular 
free density of states p°(e) = — {%) 2 - The half-bandwidth D corresponds to a hopping parameter t between 

neighboring d— sites with D = 2nli. We set D equal to unity. This density of states has the further advantage that 
naturally incorporates the same band-edge behavior as a 3— dimensional model. For a discussion of the hypercubic 
lattice see Ref. 12 



As a result, we arrive at the following local effective action 



S local =-E/ I dTdT ' ^Ar)g^(r-r')Mr') + dr (n /T (r) - ~)(n n (r) - ~) (5) 

where ip a = {f a , d a } is the operator associated to a site that we call origin (of a given sublattice, say, A). The local 
inverse propagator Q^ 1 has the explicit form 

i f iui V \ 

^WM= T/ ,, , D 2 <n i (6) 



V iu- ^-[g ctB WM 

where G is the "cavity" Green function which has the information of the response of the lattice. In the present case 
of a Bethe lattice it simply becomes G — G. The symmetry properties of the two sublattice representation imply 

G a A = G- a B (7) 

so the self-consistency condition can be more concisely rewritten as 
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D 2 

[0 o <taW(H =iu- —[G-<y A ]dd(iuj) (8) 
and the lattice index can be dropped. 

At this point it is worth pointing out that Go should not be confused with Go- While the latter is the free 
propagator of the lattice model, the former corresponds to the free propagator of the effective impurity model which 
is local. Furthermore, Go is a quantity that is unknown a priori, and that has to be solved for self-consistently from 
equations (§|,|). i— jj— ■ 

We solve this system of equations numerically by a recently introduced exact diagonalization procedure.li3~li3 It 
basically consist in introducing a parametrization for the Gdd i n terms of the parameters of two continued fraction 
expansions that describe the "particle" (ui > 0) and "hole" (u < 0) excitations. 

G drf M = G>H + G<>), (9) 

with 

(gs\dtf\gs) 



uj — a. 



o 



<nM - — fr'i^g, — do) 

lo - a - — 



u— a, — 



1 ,,-„<- 



The continued fractions can be then thought of as resulting from an effective electronic bath where the impurity site 
is embedded. This bath consists of two chains, one for each continued fraction. The hopping elements between sites 
of the chains are the parameters bf and the atomic energies of the sites are the af^. In Fig. |l| we schematically 
show the resulting effective hamiltonian model for the local impurity site plus the effective electron bath. It can be 
straightforwardly verified that it renders the self-consistency equation (ra) . 

Swithching on the interaction U, the cluster of Ns sites is then exactly diagonalized and the local Green functions 
are calculated as continued fractions. The parameters of this newly obtained continued fractions, are then feed back 
as new parameters for the effective bath chains. It should be noted that in the present case, where solutions with 
antiferromagnetic order are allowed, one should connect the spin a chains to the spin — er local d— site as follows from 
equation (||). 

The self-consistency is thus translated into the self-consistent determination of the parameters of a continued 
fraction representation of Gdd- This numerical procedure represents an essentially exact solution of the model. The 
only approximation consists in the truncation of the length of the continued fractions, due to the finite size of the 
effective electron bath that can be dealt with. It has been demonstrated elsewhere that this approach is in excellent 
agreement with the solution obtained from a quantum Monte Carlo calculationO Moreover, it allows the investigation 
of regions of parameter space that are inaccessible for QMC, for instance, the T — > and large U limit. 



III. MAGNETIC ORDER 



A. Phase Diagram 



We have performed extensive calculations for effective hamiltonian clusters of Ns = 4, 6, and 8 sites in the (U, V) 
parameter plane for fix D = 1. In Fig. || we show the exact magnetic phase diagram of the model. The parameter 
space is divided into two regions, one paramagnetic and the other with long range antiferromagnetic order. As is 
shown in the inset, the boundary line has the approximate form U c oc V 2 . The interpretation of this result is simple. 
The onset of magnetic order is controlled by the crossing of two energy scales that depend on a single parameter J 
proportional to the exchange constant If >> 1, the magnetic moment on the /—site is Kondo-quenched by 
the band of d— electrons, rendering a paramagnetic state. On the other hand, for -j^ << 1, the RKKY interaction 

(~ (-j) Q , a > 1) becomes dominant over the Kondo effect (~ e~^^). Therefore, the screening of the local moment 
becomes incomplete and the residual magnetic interaction along with the bipartite nature of the lattice drives the 
system to an antiferromagnetic state. We thus see that these energy scales cross at a critical value J c , and in 
consequence the (U, V) parameter space is split into two phases with a boundary that obeys ^j- J c . We find from 
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the numerical solution that J c ps 0.0751?. This result-is consistent with Doniach's estimate for the critical value of 
J c . In the context of a one dimensional Kondo latticeO, he found that at the mean field level jj ~ 0(1) and argued 
that fluctuation effects should strongly reduce the ratio. 

It is interesting to note that this result is qualitatively.-similar to the recently obtained by Vekic et al. for a 2- 
dimensional finite size lattice using quantum Monte CarloH This lack of dimensional dependence represents a strong 
evidence for the applicability of the LISA method in the study of electronic systems where localized states and 
correlations play an essential role. 

As it turns out, the transition is of 2 nd order along the whole critical line. This will be illustrated by the study of 
the magnetization below. 

B. Staggered Magnetization 

We now consider the evaluation of the staggered magnetization. Since the present effective cluster hamiltonian 
method treats explicitly the / and d— sites, their local magnetization can be simply and efficiently evaluated as 
expectation values on the ground state. 

m Z f = {gs\n s] - n fl \gs), m Z d = [gs\n d ] - n dl \gs) (11) 

In Fig. |3|, we show the magnetic moments as a function of the hybridization V for different values of the interaction 
U. Both moments mzd and mzf vanish continuously at the transition for any value of U indicating that the transition 
is of 2 nd order along the whole critical line. 

A notable feature of the results is that the size of the magnetic moments of the / and d— sites are very different, 
except close to the critical point. We also point out that while the magnetization of the /—site becomes always rapidly 
saturated at mzf ~ 1, independently of the strength of the interaction, the same does not occur on the d— sites. As 
U is increased, the maximum of mzd is enhanced. This suggests that it would be interesting to consider a limiting 
case of the model with both V, U — ► oo (or equivalently D — ► 0) and ^j- kept constant (i.e., — * oo). This model 
may present an almost compensated total magnetic moment. 

An imp|Oj|tant remark is that the deep Kondo regime, with an exponentially small gap in the paramagnetic density 
of statespB'ta is actually hindered by the onset of the magnetic instability. This is relevant for the understanding of 
the very small gaps observed in the Kondo insulators. This situatipa could, in principle, be modified by introducing 
frustration in the model, as for instance through n.n.n. hoppingEll, or by the consideration of degeneracy in the 
/—orbital. 

Let us now consider in more detail the behavior of these quantities in the case where the critical line is approached 
and also when J — > ► 0). 

In Fig. H we show the magnetic moments as a function of the hybridization V for a fixed value of the interaction 
U = 2 in clusters of various sizes. A similar result was obtained by Sun et al, using a variational Gutzwiller 
approximation.!^ 

As was already pointed out in that paper, the behavior of the magnetization as V — > and as V — > V c , near the 
phase transition, is quite different. In the latter case, both mzf and mzd vanish when the Kondo effect becomes 
relevant close to the transition line. In the former case, however, we see that while the magnetic moment of the /—site 
becomes saturated and fully developed rrizj = 1, the local moment of the d— electron band goes to zero. This result 
is demonstrated by the data shown in the inset, where we plot the magnetization at at small V = 0.01 as a function 
of the size of the cluster. 

It is also important to note from the comparison of the results from clusters of size N$ = 4, 6 and 8 shown in 
the main part of the figure, that although the value of the magnetization is notably affected by the system size, 
the the critical point for the ordering transition remains basically size independent. This feature makes the phase 

diagram presented in the previous section an essentially exact result for this model. Also, a scaling oc ^Jl — p^for 

the moments mzd and mzf is found near the critical lineE! Moreover, we note that the results for mzf, which are 
essentially independent of the system size, can be very well fit in the whole interval < V < V c with the function 

■\Jl — with only V c as fitting parameter. 

We now investigate the behavior of the moment along the other direction in the (U, V) plane. In Fig. |5| we show 
the results for the magnetization as a function of U with fixed V — 0.5. A notable feature is that a relatively large 
magnetization on the d— sites seems to persist for extremely high values of the interaction U ~ 100 (which are easily 
accessible within the present numerical technique), despite the fact that the RKKY coupling should be very small. 
However, as the inset demonstrates, when the data are plotted against the relevant energy scale of the model, i.e. 
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J ~ ^j-, the behavior turns out to be qualitatively similar as in the previous case (which is also an indication that 
finite size effects are important in this parameter region). 

To fully elucidate the behavior of the magnetization in the limit ^j- << D < U (J/D — » 0), we obtain the Hartree- 
Fock solution of the model (c.f. Appendix). The approximation becomes essentially an exact result in this limit, 
where correlation effects arc less important. We find that the staggered magnetization of the d— sites behaves as 

m ^^ Zn (^r)> J<<D ( 12 ) 

V/ 2 

where J = — . This result is plotted for comparison with the corresponding data from exact diagonalization in 
Fig. | (dotted line). 

To end this section we would like to briefly consider the behavior of the double occupation (D). This quantity is 
simply related to the square of the magnetization by (m 2 ) = 1 — 2(D). The results for fixed V and as a function of U 
from exact diagonalization are presented in Fig. [| For comparison, we also plot similar results for the two-site local 
model which corresponds to the limit where the band- width D — ► and can be solved analytically. When U = 0, 
we have (D) = 1/4 for both the d and /—sites since, in the non-interacting case, the four spin configurations at 
each site appear with equal probability. As U is increased, the singlet state becomes energetically favorable and the 
double occupation decreases. This behavior is qualitively well captured by the two-site model, however, as it does 
not contain the lattice feedback, it underestimates the value of (Df) and overestimates (Dd). This difference can be 
understood from the following arguments: For the /—electrons the lattice contains the RKKY interaction (that the 
two-site model lacks) which favors the development of the local magnetic moment to minimize the energy. On the 
other hand, for the d— electrons, the enhancement of the double occupation respect to the two-site model is simply 
due to their band-like character which is brought in through the self-consistency condition. 

When the interaction U is further increased the transition into the antifcrromagnctic state takes place, as shown 
by the results for the staggered magnetization mzd that we reproduce in the figure for comparison. The most notable 
feature is that (Dd) has non-monotonic behavior. It experiences an upturn and returns to the non-interacting limit 
value when U — ► oo. This behavior merely reflects the fact that as U increases, the effective magnetic exchange 
interaction J between the sites becomes small. More precisely, one may perform a Schrieffer- Wolff transformation! 2 .!] 
and realize that d and /—sites become gradually decoupled with the d— electrons approaching the non-interacting 

band limit, as J ~ ► 0, i.e., when U — ► oo. On the other hand, for the case of the /—electrons, where the local 

repulsion is the dominant interaction, the simple two-site model can still capture the correct qualitative behavior. 



IV. ITINERANT AND LOCALIZED MAGNETIC BEHAVIOR 



In this section we consider the behavior of the /—electrons as they evolv&irom the itinerant to the localized magnetic 
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regime. Different aspects of this question were first discussed by Doniachli-3 and EvansEj, and later by Sun et 

The existence of two qualitatively different regimes was already suggested by the two routes in which the system 
looses its magnetic order: (i) when V — > or U — > oo (J — ► 0), the /—moments became free local spins with mzf —> 1 
and the d— electrons, having mzd — * 0, approach a non-interacting band; and (ii) at the 2 nd order critical line, as 
both mzd and mzf — > 0. 

In what follows, we shall gain new insights on these two distinct regimes by studying the redistribution of spectral 
weight of the electronic density of states. 

At T = 0, the Green functions of the problem can be obtained as two continued fraction G ca (uj) — G^^ + Gf^^), 
with G> CT and G< a corresponding to "particle" and "hole" excitations respectivelyEZrEj 



Gca(z) = {9s\Ca TTZ ^rC^gs) 

Z - (tieff ~ tigs) 

G<{z) - (gs\ct 1 — c a \gs) (13) 

Z + ( H eff - ti gs ) 

where the anihilation (creation) operator c (c^) generically stands for d {df) and / (/^), and H e ff is the effective 
cluster hamiltonian introduced before. The —a Green function is by symmetry G c - a (uj) = — [G ca (— ui)\* . 

In Fig. [7] we display the density of states for the / and d— electrons. They are obtained from the imaginary part of 
the local Green functions of clusters with 8 sites. Note that, unlike the quantum Monte Carlo technique that has also 
been applied to this problem, our results are directly obtained on the real axis with no need for analytic continuation. 
However, the price we pay is a discrete number of poles that we broaden in the figures for easier visualization. 
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The figure illustrates the changes that take place in the density of states as we move through the different regimes 
on the phase diagram. The three sets of curves correspond to the antiferromagnetic phase, as we move away from the 
2 nd order critical line towards the small J regime (V 2 << UD). The two curves on the top of the figure correspond 
to U = 3 and V = 0.5 that places the system close to the critical point. Both spectra are slightly asymmetric which 
correspond to the simultaneous vanishing of the magnetic moments that we discussed before. A notable aspect is that 
a central feature can be clearly observed in the /—electron spectral function. This corresponds to the Kondo coherent 
quasiparticle excitations that are characteristic in the paramagnetic solution of the model. In this case it is interesting 
to observe that the quasiparticle excitations can survive the onset of magnetic order. In Fig. || we display the behavior 
of the /—electron density of states in greater detail. The symmetric curve (grey line) is obtained in the paramagnetic 
region very close to phase boundary. This solution corresponds to the well known scenario of a hybridization band 
insulating state with the hybridization amplitude being renormalized by the effect of the correlations. The low energy 
features are originated in the splitting of the Kondo resonance due to the lattice effect. Sjhis scenario is borne out 
from many-different approaches, as for instance the Gutzwiller variational wavefunctionHu, large iVEm, and large 
dimensionst3. We now proceed to slowly move into the antiferromagnetic phase. The comparison of the paramagnetic 
solution to the curves obtained close to phase boundary, but in the antiferromagnetic regime, provides us new insights 
as we realize that the transfer of weight actually occurs at all energy scales. Both, the high energy features and the 
low energy excitations associated with the quasiparticles on the left side of the spectrum, are loosing their weight 
at a similar rate. This interesting behavior actually makes difficult the task of reducing the the problem to simpler 
effective model that could be more tractable. It also does not allow for a simple interpretation of which aspect, the 
itinerant quasiparticles or the localized moments, is the driving force behind the transition. 

We now come back to Fig. as we decrease the coupling between sites, setting U = 3 and V = 0.3 (central curves). 
We, thus, move away from the 2 nd order line well into the antiferromagnetic regime. The quasiparticles features have 
gradually disappeared as the weight is transferred to the high energy part of the spectrum. This corresponds to the 
rapid increase in the size of the /—magnetic moment as V 2 becomes much greater than UD and the /—electrons 
localize. 

As we finally enter the small J regime, setting V = 0.1 (lower curves), we clearly observe how the electron become 
decoupled. The /—electron spectrum corresponds to a fully saturated and polarized moment with the totality of the 
spectral weight at high energies. The lack of any low frequency feature indicates that these electron have essentially 
localized. On the other hand, the d— electron density of states seems to have closely approached that of a non- 
interacting band. However, we shall see next, this is not exactly the case. 

Once more we can gain further insights in this limit by using the Hartree-Fock approximation, which is very accurate 
as J << D (V 2 « UD). We focus on the d— electrons which have a more interesting behavior. In this case, we can 
obtain an analytic expression for the density of states (c./. Appendix), 
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We plot this density of states, along with the results from an 8 site cluster, is Fig. |9|. It is remarkable to observe 
the effort that the few poles from the exact diagonalization method do to best approximate the essentially exact 
continuum result. Even details as the small weight at w ~ 1.5 that originates from the hybridization with the high 
energy feature of pf (c.f. Fig. [7]) are well captured. As expected, the d— electron density of states is insulating with 
a small gap of size of the order of J. However, it is surprisingly far from regularly approaching the non-interacting 
limit. A very sharp peak develops on the low frequency edge of the gap. When the J — > limit is approached 
the peak becomes narrower and gradually looses its spectral weight. This feature could, in principle, be observed in 
infrared photoemission spectroscopy. Also, in transport measurements it may be possible to observe a strong pressure 
dependence at temperatures of the size of the gap. 

V. CONCLUSIONS 

In this work we have considered solutions with magnetic order in the periodic Anderson model in the limit of 
large dimensionality using the LISA method. Our results validate and expand on that of Sun et al. who obtained 
approximate solutions based on a variational procedure. 

The T = exact magnetic phase diagram of the model is obtained in detail by a recently introduced numerical 
algorithm that is based on the exact diagonalization of an effective cluster hamiltonian. We find a paramagnetic- 
antiferromagnetic 2 nd order phase boundary that splits the (U, V) parameter space in two regions. The critical line 
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obeys U c oc V c 2 as U becomes large. This is a concrete realization of earlier ideas of Doniach in the context of the 
Kondo lattice, and is consistent with recent numerical results on 2-dimensional finite lattices. 

The consistent picture that emerges from this work, very similar to the two dimensional lattice results, represents 
further evidence of the applicability of the LISA method to investigate the physical behavior of systems that are 
characterized by the simultaneous presence of localized orbitals and strong correlations. 

From the determination of the phase diagram we found that the paramagnetic Kondo regime with an exponentially 
small gap, which is relevant for the description of Kondo insulator compounds, is inaccessible in all parameter space 
due to the onset of the antiferromagnetic instability. 

The calculation of the magnetic moments in the phase with long range order, revealed that the possibility of a 
nearly perfect magnetic spin compensation that would be relevant for the problem of the "small moments" , is not a 
general feature of the solutions of the model. 

An almost perfect compensation of the moments is only obtained near the 2 nd order boundary line, where the 
total site magnetization starts to become non-zero as the Kondo screening ceases to be complete. Therefore, a small 
moment regime would require a fine tunning of the parameters of the model, which is not a priori a feature that one 
expects to find in real systems. The question weather this situation would persist in the non particlc-holc symmetric 
case and how the introduction of frustration may modify the present picture is currently under investigation. 

Finally, the method allows to obtain reliable information on the distribution of the spectral weight in the density 
of states which provided interesting new insights on the nature of the phase transition and that of the different 
antiferromagnetic regimes. 

ACKNOWLEDGMENTS 

Useful discussions with G. Kotliar, E. Miranda, A. Georges and L. Laloux are gratefully acknowledged. 



APPENDIX A: THE HARTREE-FOCK SOLUTION 

We present in this appendix the solution of the self-consistent equations (||,^],||) in the limit ^j- << D < U. It is 
a useful exercise that would allow the interested reader to better grasp some details of the method. Performing the 
Hartree-Fock decoupling to interaction term of the effective local action (||), and using that {rif a — nf_ a ) — (mzf) 
and {nf a -\-nf_ a ) = 1, the action becomes quadratic and we can immediately obtain the local Green function for the 
J (j electrons . 

G Af± „ M = — jj \ pr (Al) 

where A and B denote the sublattice. We now combine this expression with the symmetry relation between the 
sublattices Green functions (0), to obtain 

G Af± „ (w) - — v 1 — (A2) 

to =F mzf y - V 2 G A d^ a (w) 

since the expression involves only quantities on the same sublattice, the index may be dropped. We now need to 
find an expression for the d— electron Green function. Once more, from the local action we obtain, 

G Ad±a M = y2 1 — (A3) 

where from the same symmetry considerations we have, 

G d± „ H = v , \ 2r . : (A4) 

This equation is actually two coupled equations for Gd T and . We then replace Gd„ into its same expression, in 
order to decouple them. We obtain, 
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Gd ± A")= V 2 ~ 2 i (A5) 

We are left now with two quadratic equations, each with a single unknown, which we can solve to obtain the explicit 
result, 



G d „ M = 1 - 2 — ( AO ) 

from which follows the d— electron density of states that we express in terms of the half-bandwidth D = 2t 



Pd a M 



irD 2 



\ 



v 2 



v 2 



-D 2 



V 2 

UJ - mzf y 



V 2 



We still need an equation for mzf. It is simply obtained from its definition, 



(A7) 



1 

mzf = — 



Im(G n (ui) - G fl (ui)) du 



(A8) 



Thus, (A2,A^.A8) form a system of equations that have to be self-consistently solved for. 

In general, these equations can be solved numerically. However, from the results of section III we learn that in the 
limit ^j- << D < U, which we are currently concerned with, we can take mzf = 1 to an excellent approximation (c.f. 
Fig. |). Thus, equation (A7) immediately leads to the expression for the density of states ( |i"4| ) in section IV. We can 
finally obtain the behavior of the magnetization mzd by replacing into 



mzd 



Im(G d ^(uj) - G dl (uj)) da; 



(A9) 



an extracting the leading contribution in this limit (c.f. Eq. (h_2) in section III). 
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FIG. 1. Effective cluster hamiltonian. It consists of a central impurity site with a d and /—electron site connected to an 
effective electron bath. U is the on site repulsion of the /—site, V is the d — f hybridization matrix element, and t is the d — d 
hopping amplitude. 



FIG. 2. Phase diagram of the PAM. The inset shows in a log — log plot that the 2 nd order critical line obeys U c ~ V c 2 as U 
becomes large. Ns = 6. 

FIG. 3. The staggered magnetic moments m Z d and m Z f as a function of the hybridization V for different values of the 
interaction U = 1, 5, 10. Ns = 6. 



FIG. 4. The staggered magnetic moments mzd and mzf a function of the hybridization V for a fix value of the interaction 
(7 = 2. For comparison we plot similar results for Ns = 6 (bold lines) and Ns = 4 (thin lines). Note that while the magnitude 
of mz depends on the size of the cluster, the position of the critical point does not. The dotted line in the top part of the 



figure indicates the Hartree-Fock results, while the one in the bottom is the fit y 1 — {y^) ■ The inset shows mzd and mzf 
(top and bottom) for a small value of V = 0.01 as a function of the inverse of the cluster size l/iVg. m Z d becomes small while 
mzf remains saturated at mzf = 1. 

FIG. 5. The staggered magnetic moments mzd and mzf as a function of the interaction U for a fixed value of the 
hybridization V = 0.5. Ns = 6. The inset contains the same quantities plotted as a function of V (for fixed U = 5), and of 
1/(7 (for fixed V = 0.5). As the relevant energy scale of the problem is J ~ V 2 /U, the behavior of the magnetic moments is 
qualitatively similar in both cases. 

FIG. 6. Double occupation (D) as a function of U for V = 0.5 from exact diagonalization of a 6 sites cluster (bold lines). 
(Dd) is the upper curve and (Df) the lower one. For comparison we include in dotted line similar results for a 2 sites cluster 
((Dd) = (Df)). The thin line corresponds to m Z d that indicates the magnetic regimes. 

FIG. 7. The local density of states P/„,Pd„ (left and right) in different magnetic regimes. The top spectra correspond to 
V = 0.5 and U = 3 close to the 2 nd order critical line. The low energy features of the P<j f correspond to quasiparticle AFM. 
The bottom spectra correspond to V = 0.1 and [7 = 3. The single peak at a frequency ~ y and the absence of quasiparticles 
indicate local moment magnetism. The middle spectra correspond to V = 0.3 and U = 3 in the crossover region. Ns — 8. 

FIG. 8. The local density of states pf a for U = 2 and V = 0.45 in the paramagnetic phase (grey line) and V = 0.40, 0.35 
in the antiferromagnetic phase (dotted and full lines). 

FIG. 9. The local density of states pd a for (7 = 3 and V = 0.2. The bold line corresponds to the Hartree-Fock solution and 
the thin line to the exact diagonalization of an 8 sites cluster. 
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